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1 Introduction 



The generalization of the Gell-Mann-Low or adiabatic theorem [|]] (see also 0) was es- 
tablished recently pi by one of the authors, with the intention of making a certain class 
of "non-perturbative" phenomena accessible to perturbative methods. The emergence of 
bound states is a perfect paradigm for such phenomena: from a Lagrangian field theory 
point of view, bound states are characterized by poles in the S-matrix and as such can only 
be generated by summing an infinite set of Feynman diagrams ||. The most traditional 
method is to choose a particular (finite) set of diagrams contributing to the four-point func- 
tion, and to iterate these diagrams to all orders via an integral equation, as, for example, 
the Bethe-Salpeter equation JJ, [J (for a review see Ref. 0). In that case all two-particle 
reducible products of the initial set are included in the final result. Therefore the initial set 
itself should contain only two-particle irreducible diagrams. In most calculations only the 
one-particle exchange diagram is taken into account in the kernel of the integral equation. 
The iteration of this kernel through the solution of the integral equation leads to the set of 
so-called ladder diagrams, and the corresponding approximation to the four-point function 
is known as the ladder approximation. 

The Bethe-Salpeter equation, or rather the Bethe-Salpeter approach, has many problems 
0, ||, most of which are related to the appearance of a relative-time coordinate. The 
covariant, Lagrangian, or space-time formulation of field theory was originally designed 
for scattering problems, where the interaction distances are short and consequently the 
interaction time is short as well. However, for a bound state, the "interaction time" is 
infinite, and for massless exchange particles the interaction range is large. For such problems 
the Hamiltonian, time-independent approach is better suited, but difficult to implement 
consistently. 

Since the formulation of the Bethe-Salpeter equation, there has been a wealth of alterna- 
tive proposals for bound-state equations in quantum field theory. In the beginning, despite 
the drawbacks mentioned above, they all started from the Bethe-Salpeter equation, with 
the relative time eliminated in one way or another. These equations are generically known 
as quasipotential equations 0, with the Blankenbecler-Sugar equation || || and the Gross 
or spectator equation jllj as well-known examples. However, even though these equations 
have to satisfy some unitarity conditions, much arbitrariness is left, and results depend on 
the way the four-dimensional equations are reduced to three-dimensional ones [TT| . 

It is important to notice that all these approaches are rooted in the Bethe-Salpeter equa- 
tion and therefore in Lagrangian perturbation theory. For the description of physical states, 
and in particular bound states, as mentioned before, the Hamiltonian formalism is prefer- 
able over the Lagrangian one, the tool of choice being the time-independent Schrodinger 
equation. Indeed, it is clear that in a Hamiltonian approach a relative-time coordinate has 
no place, thus eliminating one of the major problems of the Lagrangian approach right from 
the start. 

More recently, other routes to bound-state equations have been explored. Here we just 
mention a few. Firstly, there is the Feynman-Schwinger Representation (FSR) approach 
T2|]-|14j], which starts from the path-integral representation of the four-point function. The 



integral is performed partly analytically and partly with Monte Carlo techniques in Euclidean 
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space, hence the FSR is closely related to the Lagrangian approach. However, most problems 
associated with the Bethe-Salpeter equation are absent, or not yet discovered, as it is hard 
to recover excited states in this approach. 



Secondly, there are a number of light-front field theory approaches p5|-p3[, which either 
work in discrete momentum space, or at low orders in a Fock-space expansion, but cover 
a wide range of approaches from rather traditional to innovative. In all these Hamiltonian 
approaches, (approximations to) physical states are calculated. 

Thirdly, a few investigations have been made using the Haag expansion, or N-quantum 
approach P^j] , but, probably because of computational difficulties, it has never been suc- 
cessful. Fourthly, one of the authors has recently explored the same problem in an ordinary 



equal-time Hamiltonian formulation with a Fock-space truncation ||25|| . Self-energy effects 
have been included, which had been ignored in earlier work. Finally, the other author has 
applied Regge theory to extract bound state energies from the leading Regge trajectories 
as calculated by one-loop renormalization group improved Lagrangian perturbation theory 



However, few of these approaches formulate a perturbative approximation to an effective 
Hamiltonian, which can then form the basis of a bound-state equation in a given particle 
sector. This and the beforegoing considerations have been the major motivation for the gen- 
eralization of the Gell-Mann-Low theorem, which may be considered to be a particularly ef- 
ficient formulation of Hamiltonian perturbation theory. Other non-perturbative phenomena 
like spontaneous symmetry breaking and vacuum condensation are hoped to be describable 
by the same method. In each case, the application of the Gell-Mann-Low theorem leads to 
a series of perturbative approximations to a certain equation. In analogy with the integral 
equations of the traditional Lagrangian approach, it is the solution of this equation (or any 
of its approximations) which leads to the wanted non-perturbative information. In the case 
of the bound-state problem, this equation is nothing but the time-independent Schrodinger 
equation, where the potential part is given in the form of a perturbative series. 

Unfortunately, in this Hamiltonian approach there is also a counterpart to the arbitrari- 
ness in the Lagrangian approach to bound states: different mappings to the unperturbed 
subspace that appears in the generalization of the Gell-Mann-Low theorem, lead to non- 
equivalent Schrodinger equations at finite orders of the perturbative expansion. A future 
discussion of this issue will take properties like renormalizability, and Lorentz and gauge 
invariance at finite orders of the perturbative expansion into account, but it will also have 
to consider applications to simple models as in the present contribution. For the time be- 
ing, we will focus on the particular choice discussed in Ref. @, which has somewhat unique 
properties. 

In this paper we will demonstrate the applicability of the generalized Gell-Mann-Low 
theorem to bound-state problems with the help of a very simple model, and compare with 
the Bethe-Salpeter and related approaches, arguing the superiority of the former method 
over the latter. The problem of choice is the calculation of bound states of two distinguishable 
scalar constituents interacting through the exchange of a third scalar particle. In the case of a 
massless exchange particle, and in the ladder approximation to the Bethe-Salpeter equation, 
this model is known as the Wick-Cutkosky model |2~7fl . 

In the following section, we will derive the effective Schrodinger equation for this scalar 
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model from the generalized Gell-Mann-Low theorem to lowest non-trivial order, emphasizing 
the diagrammatic representation of the effective Hamiltonian. The resulting diagrams are 
similar to Goldstone diagrams but, unlike the latter, in general they cannot be combined 
into on-shell Feynman diagrams. The third section presents an analytical study of the 
non-relativistic and one-body limits, demonstrating the consistency of the formalism in 
these particular cases. In the fourth section, we solve the effective Schrodinger equation 
numerically, establishing that all solutions are physical and confirming the analytical results 
of the preceding section. Finally, the last section presents a critical discussion of the results 
and a comparison with the Bethe-Salpeter equation and other approaches. 

2 The Schrodinger Equation 

The new method we propose and explore here for the solution of the bound-state problem 
in quantum field theory, is a generalization of the Gell-Mann-Low theorem || . The essential 
statement of the generalized Gell-Mann-Low theorem (GGL for short) is that the adiabatic 
evolution operator maps linear subspaces invariant under the free Hamiltonian H to linear 
subspaces that are invariant under the full Hamiltonian H = H Q + Hj. The full Hamiltonian 
is then diagonalizable in that subspace and one need not consider the eigenvalue problem in 
the whole Fock space. The mapping between the two linear subspaces considered in Refs. 
|| is such that when followed by the orthogonal projection from the ii-invariant subspace 
back to the Ho-invahant subspace, it gives the identical map in the latter. In order to assure 
this normalization, the adiabatic evolution operator has to be accompanied by a mapping 
within the if -invariant subspace, and the combination of these two maps has been called 
the Bloch- Wilson operator. 

The proof presented in Ref. rests on the existence of a perturbative series for the Bloch- 
Wilson operator, or more precisely on the existence of its adiabatic limit. It should not be 
taken for granted, but rather depends on a sensible choice for the if ~i nvar i an t subspace. 
In practice, one will only calculate to a certain perturbative order, and the existence of the 
limit is then only guaranteed (or needed) up to that order. The explicit calculation to be 
performed in the present work will demonstrate that a reasonable approximation can be 
achieved by this procedure to lowest order, and that indeed a sensible and natural choice 
(to this order) for the if ~i nvar i an t subspace exists. It should also be noted that the Bloch- 
Wilson operator establishes a similarity transformation between the two linear subspaces 
so that the problem of diagonalizing the Hamiltonian H in the if -invariant subspace can 
be transported back to the i?o~invariant subspace. The image of H under the similarity 
transformation in that subspace is Hbw, the Bloch- Wilson Hamiltonian. Below, Hbw will 
be calculated to lowest non-trivial order for the problem at hand. 

We will now introduce our notations and also restate the basic formulas (for their deriva- 
tion and discussion see Ref. ||). The adiabatic evolution operator U e is given by the series 
expansion 
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dt n e-<^+- + ^T[H I (t 1 )---H I (t n )\, 



(1) 



3 



where 



Hj(t) 



JH t 



H ie 



-iH t 



(2) 



The Bloch- Wilson operator U bw is defined by the adiabatic limit 



U BW = lim U e (P U t P 



,-1 



(3) 



Pq being the orthogonal projection to the if -invariant subspace Qo- Granted the existence 
of Ubwi the Bloch-Wilson Hamiltonian is simply 



Hbw — PqHUbw ■ 



(4) 



The scalar model we consider in the following consists of three scalar particles, two of 
them carrying some charge and one neutral, with masses tua, n^B and fi, respectively. The 
two charged particles are coupled to the neutral one via a charge-conserving cubic coupling 
with strength g. The coupling constant g has the dimension of a mass, and for our purposes 
we might as well have chosen two different constants g^ and gs for the coupling to the two 
charged fields. The Hamiltonian of the model is 



H 



H 



Ho + Hj 



d x 



A(x)K " V 2 )0a(x) + 4(x)(m| - V 2 )<Mx) 



+ -^(x)( / , 2 -V 2 Mx) 



Hi 



(fx 



(x)0a(x)^(x) + 



>B 



x)v?(x) 



(5) 



where we have normal-ordered the interaction part, mainly for reasons of presentation. It is 
well-known that the perturbative vacuum is unstable in this model field theory. However, 
perturbation theory around this vacuum is well-defined order by order, and no instability 
appears as long as the number of (perturbative) constituent particles in any intermediate 
state is finite. This model has been much used in the past to test bound-state equations. 
In particular, it has played an important role for the Bethe-Salpeter equation, where the 
ladder approximation in the special case \i = is known as the Wick-Cutkosky model ||27| . 
The bound states usually considered are the ones containing one A- and one S-particle as 
constituents. The model is then the simplest one can possibly think of. 

The Fock space will be built up as usual by repeated application of the creation operators 
to the vacuum {Q\Q) = 1. We will normalize the creation and annihilation operators in 
momentum space in such a way that 



[a l (p),al(p')] = (2 7 r) 3 5(p-p' 



(6) 



for all three particles (labelled by the index i). This normalization is of course not relativis- 
tically invariant, but it turns out to be convenient, and is closely related to the quantum- 
mechanical wave functions. Besides, Lorentz invariance is not manifest in a Hamiltonian 
approach anyway. 
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We can now define the if ~i nvar i an t subspace Sl we wn l use for the application of the 
GGL. Considering that we want to determine bound states of an A-and a P-particle which 
will contain components with different momenta of these particles, a natural choice is 

tt = span jlpAjPs) = a] i (p A )a t B (p B )\^) P^PbSK 3 } , (7) 
and the corresponding orthogonal projector reads 

f d 3 p A d 3 p B 

Po = 7 (2^(2^ |PA ' Ps)(PA ' Ps| - (8) 

We could have restricted Qq to states with fixed total momentum P = p^ + p#, but it is 
more instructive to see how this constraint arises from momentum conservation. 

The rest of this section will consist in calculating and interpreting the Bloch-Wilson 
Hamiltonian Hbw to lowest non-trivial order. In Eqs. (|3]) and (Q), Hbw was defined as the 
adiabatic limit e — > of the operator 



-i 



PoHU^PoU^Po)- 1 
= PoHoU^PoU^)- 1 + P Q H I U e (P U e P ] 
= HoPo + PoHjU^PoU^o)- 1 , (9) 

where in the last step we have made use of the fact that Qq is a direct sum of eigenspaces 
of H , hence PqHq = H P Q , and of the fact that one can insert an operator P between 
U e and (PoC/ e Po) _1 . To first order in Hj, one needs to keep only the trivial term for U e , 
i.e. U e — 1 + O(Hj). But since Hj necessarily changes the number of particles, one has 
P HjP = 0, and Hbw reduces to the non-interacting part, Hbw = H P . 

Consequently, for the lowest non-trivial order one has to expand U e to first order in Hj, 



U e = l-i / dte-^H^t) + 0{Hj) 

J —oo 



(10) 



Using PqHjP = as before, the Bloch-Wilson Hamiltonian to lowest non-trivial order is 
given by 



H BW = H P - % \ dt e'^PoHj (0)i?,(f) P + O(Hj) , (11) 

J —oo 

where the limit e — > is understood. 

We could now go on and directly determine an explicit expression for Hbw to this order, 
but it is conceptually clearer (and in principle necessary) to consider the zero- and one- 
particle sectors first. A posteriori we will see how an analysis of the diagrams contributing 
to the Hbw of Eq. (|TT|) yields the same information. To begin with, we study the vacuum 
sector to the same order in the coupling constant g. The Hamiltonian in Eq. ( pf) then 
remains unchanged except for the replacement of Pq with the vacuum projector 

P n = \Q)(Q\. (12) 
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t = t = 

A B 



Figure 1: The second-order contributions to the vacuum energy according to Eq. fll4]) . The 
propagators of the charged A- and 5-particles are represented by thick and thin lines, 
respectively, while for the neutral particle we have drawn broken lines. Integration over t 
from — oo to for the vertex to the left and integration over the spatial position x for both 
vertices is understood. 



We can simplify Eq. (O) to 

H BW = (e -i ^ dte-MiniTtfjWHjWlfy) P n + 0(Hf), (13) 

where Eq = (Q\Hq\Q) is the vacuum energy, due to the fact that we have not normal- 
ordered H Q (there would be an additional contribution to second order in g, had we not 
normal-ordered Hj). By use of Wick's theorem, the matrix element in the second-order 
contribution to Hbw can be written in terms of Feynman propagators in position space as 



g 2 J d 3 x d 3 x'Aj(0 - t, x - x') A$(t - 0, x' - x) A F (0 - t, x - x) 
+ g 2 J d 3 x dVAf (0 - t, x - x') Af (t - 0, x' - x)A F (0 - t, x - x) 



(14) 



where 



with 



/d^k Q—ikot+ik-x 
+ * (15) 

, " ^(t)e^ ( " k '- kx) + #(-i) e ^ k '- kx n , (16) 
(27r) 3 2^ k 

w k =V/?+k*, (17) 

and analogously for the A- and 5-particles. Eq. (|TJ|) is represented diagrammatically in 
Fig- E- 

Using the non-covariant representation Eq. ([16]) of the Feynman propagators and per- 
forming the integrations over t, x, and x', the second-order contribution in Eq. (|l^) is cast 
in the alternative form 



r° 

En-E = -i dte-^^TlHj^H^t)]^) 

J — oo 



-g 2 



d 3 k d 3 k! 1 



(2tt) 3 2c^ (2tt) 3 2^, 2w k+k ,(w£ + u£, + cj k+k , 
d 3 k d 3 k! 1 
(2tt) 3 2c^ (2tt) 3 2^ 2w k+kf (u* + ug + cj k+k 
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Of course, this expression is not very meaningful as it stands. Some regularization procedure 
would have to be employed to render it finite. On the other hand, what will be of interest 
is only the difference in energy between the vacuum and other states. Since the expression 
Eq. (|i"8|) always appears as the result of evaluating the diagrams of Fig. [I], it will always 
be easy to identify and subtract it. The same expression results in usual time-independent 
Schrodinger perturbation theory and also in covariant Lagrangian perturbation theory to 
second order. 

Having treated the vacuum sector in some detail, we now go on to discuss the one- 
particle states, taking as an example the momentum eigenstates of one A-particle. The 
natural choice for the corresponding i7 -invariant subspace is Qa, the span of all these 
states, and the corresponding orthogonal projector is 

pa =S I^ |im)(pa| - (i9) 

The corresponding Bloch-Wilson Hamiltonian to lowest non-trivial order is again given 
by Eq. (11), with P A replacing P , and the relevant matrix elements of the second-order 
term are 

(pAlTlHMHrit)]^) 

= g 2 J rf 3 xrfVA^(0-t,x-x')A^(t-0,x'-x)A F (0-t,x-x / )(27r) 3 5(p A -p^) 

+ g 2 J d 3 x dVAf (0 - t, x - x') Af (t - 0, x' - x)A F (0 - t, x - x)(2n) 3 5(p A - Pa) 

+ g 2 j d 3 x d 3 x'^*(0, x) A F (0 - t, x - x') A F (0 - t, x - x>£ (t, x') 

+ g 2 J rf 3 xc/ 3 xVp;(t,x')A^(t-0,x'-x)A F (t-0,x' -x)^(0,x) , (20) 

in terms of Feynman propagators, where the wave functions ?/>p(t,x) are given by 



e p 



V£(t,x) = (21) 



minding the normalization of the one-particle states implicit in Eq. (j^). To arrive at Eq. 
(p0|), Wick's theorem has been used, including contractions like 

eu(p) tf A (t, x) = ^*(t, x) , (f> A (t, x) a A (p) = <(t, x) . (22) 



The four terms appearing in Eq. ( pD|) are represented diagrammatically in Fig. 0. The first 
two terms in Eq. (|20"D are diagonal in the basis {|pa)}, arid a comparison with Eq. ([14]) shows 
that, together with the contribution £ , (27r) 3 5(p j 4 — p' A ) to order zero, they just reproduce the 
vacuum energy Eq calculated before. On a diagrammatic level, these contributions precisely 
correspond to the unlinked diagrams in Fig. |2|. An unlinked diagram by definition contains 
a part that is not connected (or linked) to any external line. 
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t = 



t = 



B 



p'a 



+ 



PA 



P'a 



PA 



+ 



/ n ; 

V' A > * * PA 



+ 



Pa 




PA 



Figure 2: The second-order contributions to the one-particle states representing Eq. (|20|). 
The external lines with momentum labels correspond to the wave functions Eq. (PH ), repre- 
senting the free states. 



The integration over x and x' enforces 3-momentum conservation also for the linked 
diagrams in Fig. |2] (the third and fourth contribution in Eq. fl2"0p), so that these contributions 
are diagonal in the basis {|pa)} as well. Momentum conservation has another important 
consequence in the special case of one-particle states: it implies energy conservation, i.e., 
uip A = . The product of the wave functions in Eq. ( |20| ) and consequently the complete 
matrix elements then become invariant under a translation in time t. One can use this fact 
to move the right vertex in the last diagram in Fig. |2| to time —t, keeping the left vertex 
fixed at t — 0, and as a result the two linked contributions now sum up to a proper (not 
time-ordered) Feynman diagram. The whole process is represented diagrammatically in Fig. 
1 

In mathematical terms, we replace t by — t, interchange x and x' and make use of ui^ A = 
L)^i , in order to rewrite the contribution corresponding to the last diagram in Fig. |2| as 



w 2 



dte~ eltl J AdV^;(0, x)A^(0 - t,x - x')A f (0 - t,x - x')^, (f, x') . (23) 



It is hence of the same form as the contribution of the other linked diagram, except for the 
range of integration of t. Adding up these two contributions, with the use of the covariant 
representation Eq. (|T5|) of the Feynman propagators and 

dt e-^e iu,t = — — = 2vr5(cu) (24) 

_j to + ie uj — ie 

in the limit e — > 0, leads to the result 

1 



G (2) (<,PA)(27r) 3 5(p A -p^), (25) 



2 <~ l ~ Pl ' 
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Figure 3: The translation in time for the linked diagrams of Fig. |2| as a consequence of energy 
conservation implied by momentum conservation in this case. In the last Feynman diagram, 
the vertex to the right is still fixed at t — 0, while the position of the left vertex is to be 
integrated over t from — oo to oo. 



where 

G(2) (P) = ^ I , 2 * ; * 2 _,. (26) 

J (27r) 4 k z — fi 2 + it [p — k) z — m A + 

is the second-order contribution to the usual two-point function (without the 5-function for 
energy and momentum conservation). It is very important that we can identify the sum of 
the linked diagrams with a proper Feynman diagram as it appears in covariant perturbation 
theory, because it is then possible to apply the usual renormalization procedure in order to 
make sense of the (without regularization formally infinite) mathematical expressions. 

As is well-known, G^(p) actually only depends on the invariant square p 2 , hence on m A 
in Eq. (H). Putting 



Am 2 A = G^(m A ), (27) 

Ami i s °f second order in g, and the complete matrix element of Hbw up to second order 
in the one-particle sector reads 

(PA\H BW \p' A ) = ^E n + ^ A + ^p)(2ny5(p A -p' A ) + 0(g 4 ) 

= (e u + sjM\ + (27r) 3 5(p A - pi) + 0{g A ) , (28) 

where we have defined the "renormalized mass" Ma by 

M\ = m\ + Ami • (29) 

As is clear from the definition Eq. (f27|) of Ami, this renormalization corresponds to an 
on-shell scheme. To sum up the discussion of the one-particle sector, the second-order 
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contribution to the Bloch- Wilson Hamiltonian contains the correction to the vacuum energy 
and the mass renormalization to this order. It is clear that exactly the same conclusions 
hold for the sector with one 5-particle. 

We will now return to the two-particle subspace f2 defined in Eq. (0). The diagrams 
corresponding to the matrix elements of the Bloch- Wilson Hamiltonian in this subspace are 
presented in Fig. |j. They can be naturally classified into unlinked diagrams (the first two), 
linked but disconnected diagrams (the following four diagrams) and connected diagrams (the 
last two). 

It should be clear by now that the first six contributions are diagonal in the basis 
{|Pa?P_b)}) an d that, together with the zero-order contribution, they lead to the result 

(e q + s Jm\ + V \ + V^l+Pl) (2vt) 3 5( P a - p' A )(2n) 3 5(p B - p' B ) , (30) 

with Eq and M\ (and analogously M|) as defined in Eqs. (fL8D and (p9|) , respectively. In 
general, unlinked diagrams contribute to the vacuum energy, while linked but disconnected 
diagrams are concerned with the properties of the individual particles (in the case of two- 
particle states). 

As for the last two diagrams in Fig. [|, they correspond to the mathematical expression 
(Pa,Pb\V\p' a ,p' b } 

= -ig 2 f dte'W ! d 3 xc/V^;(0,x)^(0,x)A F (0-t,x-x> p 4 ;(t,x')^(t,x') 

J —oo J 

I dte-W I rf 3 xrfV^;(0,x)^(0,x)A F (0-t,x-xO<;(t,x')^(t,x') 

J —oo J 



w 2 



(31) 



with the wave functions ipp(t, x) defined as in Eq. fl2T| ) (and analogously for the 5-particles) . 
The integrations over the vertex positions x and x' enforce conservation of the total momen- 
tum Pa + Ps = Pa + Pb- Of course, in this case there is no conservation of the individual 
momenta, and so "energy conservation" is not implied, i.e., Up A + oj v 3 b ^ u)£, + u^, in 
general. Physically, this means that states with different unperturbed energies mix, which 
was obviously to be expected in the case of bound states. As a consequence, Eq. fl3l|) is 
not invariant under time translations, and there is no way to combine the corresponding 
diagrams into a proper on-shell Feynman diagram. 

We conclude that it is the connected diagrams that determine the interaction of the con- 
stituents of a bound state (in the case of two-particle states). In principle, it was necessary 
to determine the energies of the vacuum and the one-particle states to be able to identify 
that part of the energy of the two-particle states which has to do with the binding (or 
with the kinetic energy of the particles far away from the scattering region in the case of 
scattering states). However, the diagrammatic analysis allows to identify immediately the 
contributions which are related to the interaction of the constituents. 

We can now write down the Schrodinger equation corresponding to the Hamiltonian H BW 
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Figure 4: The second-order contributions in the two-particle sector. The first six diagrams 
are essentially those of Figs. [I] and @, except that they are multiplied with further propaga- 
tors which merely introduce additional 5-functions in the external momenta. The last two 
diagrams describe the interaction of the two particles to this order and correspond to the 
expression Eq. fl3T|). 
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in the two-particle sector, 



'm 2 a + p \ + ^Ml + pl) ^(Pa, Pb) 
+ / ^|^(PA,p S |^|p^p^(p^p , B ) = (S-^)^(PA,P fl ), (32) 
with the momentum-space wave function 

^(Pa, Pb) = (Pa, Pb|V>) • ( 33 ) 

By use of the non-covariant form Eq. flip] ) of the Feynman propagators, the potential Eq. 
(|3T|) can be cast in the form 

(Pa,Pb\V\p' a ,p' b ) 

9 2 



PA PB P' A P'g 



x (2tt) 3 5(pa + pb - Pa - Pb) ■ (34) 

In the latter expression, e can be put to zero, since for /i ^ the denominators are strictly 
positive, and for \x = the singularity is integrable and does not lead to any imaginary 
contribution. The singularity for /i = is of a similar nature as the one appearing in 
the usual Coulomb potential in momentum space [2S]. It is associated with a long-range 
force due to a massless exchange particle. In the spirit of perturbative renormalization, it 
is permissible to replace the "bare" masses m A and mg in the potential Eq. fl3~4|) by their 
renormalized counterparts, given that the difference between the two is of higher order in 
the coupling constant. 

One can now make use of the overall momentum conservation in order to reduce the 
dynamics to the center-of-mass system. One then has, with p = p^ = — p#, 



M\ + p^ + X /M% + p 2 ) ^(p) + J V(p, p'Mp') = m P ) , (35) 



where 



and 



E' = E - E n (36) 



v '(p-p') / ^ lix — A+^r— -)■ (37) 
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Figure 5: The first-order contributions to the physical states in the two-particle sector. The 
correspondence to mathematical expressions is as before, except that the ingoing external 
lines have to be contracted with the solutions ^(P-AiPb) of the Schrodinger equation Eq. 
(|3"2]), and the particle content of the state can be read off from the outgoing external lines. 



In this form, the equation will be solved in Section The effective Hamiltonian on the 
left-hand side of the Schrodinger equation Eq. (|35D obviously consists of a relativistic kinetic 
term and a potential term, where the operator corresponding to the potential is non-local and 
non-Hermitian. The non-relativistic and one-body limits of the potential will be discussed 
in the following section. 

Finally, we will turn our attention to the physical states in the perturbative expansion. 
They are simply given by the application of 

r° 

U BW = 1-1 dt e' e ^Hj(t) + 0(H]) (38) 

J — oo 

to the two-particle states 

d 3 p A d 3 p B 



(2vr) 3 (2tt) 3 



4>(Pa, Pb) |pa, Pb) e , (39) 



with the wave functions *P(pa,Pb) which solve the Schrodinger equation Eq. (|32|) . The 
result is easily written down and is represented in Fig. [5[ As an example, the mathematical 
expression corresponding to the last diagram in Fig. |5| is 
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W? x) ^ (t ' x) ^ (t ' x)(27r)3j(R4 - |p ^' PB ' k) 

d 3 p' A d 3 p\ 



1 /" d 3 fc 1 1 



^-|p^p' B -k,k). (40) 



a; , 



For the model considered, the operator Ubw is unitary to first order in the coupling constant. 
As a consequence, the physical states are normalized when the normalization is calculated 
to this order, provided that the wave function ^(PAjPb) is properly normalized. However, 
this is a feature that does not carry over to higher orders, and one should bear in mind that, 
on physical grounds, the criterium for a bound state is that the relative wave function ip(p) 
be normalizable, the latter describing the probability distribution of the constituents in the 
center-of-mass system. 

The operator U bw Ubw-> which appears naturally when calculating the norm of the phys- 
ical states, also plays a role when trying to make the effective Hamiltonian hermitian by a 
similarity transformation. For future investigations, we only remark here that this operator, 
when divided by its vacuum expectation value {Q\U bw Ubw\®)i leads to well-defined expres- 
sions even if one includes the second-order term arising from substituting the expression Eq. 

for Ubw- A generalization of the optical theorem applies to this case, for which the 
adiabatic e-prescription proves essential. 



3 Non— Relativist ic and One— Body Limits 

In this section, we will carry out an analytic study of the properties of the Bloch- Wilson 
Hamiltonian to lowest non-trivial order, or the corresponding Schrodinger equation Eq. (|35|) . 
A numerical solution of Eq. ( |35| ) will be presented in the next section, which at the same 
time will confirm the results to be obtained in the following. 

There are two different limits in which a formalism for bound state calculations can easily 
be tested for consistency: the non-relativistic limit, in which Eq. (|35|) should reproduce the 
corresponding non-relativistic Schrodinger equation, and the one-body limit Mb — > oo, 
where Eq. (|35|) is expected to reduce to a (relativistic) equation for particle A in a static 
potential PH. 

We will discuss the non-relativistic limit first, which is by definition the case p 2 <C 
M A ,M B , where p = p^ = —ps is the relative momentum. More precisely, this means 
that there exists a solution ip{p) of the effective Schrodinger equation which takes on non- 
negligible values only in this region of momentum space, or equivalently, which is strongly 
suppressed for p 2 > M 2 , 

M '=Ara <M ^ MB <41) 



being the usual reduced mass. If such a solution exists, the effective potential Eq. (|37Q can 
be replaced in the Schrodinger equation Eq. fl35|) by its limiting form for p 2 «C M\, M B and 
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p' 2 «Mi,Mi, 

g 2 1 4na 
p-p 

with the dimensionless effective coupling constant 

f 

16%M A M B 



^■p') = -4^^- = > + (p- P -r (42) 



a = , (43) 



which is the analogue in this model of the fine-structure constant in QED. To arrive at the 
approximation Eq. ([42]), we have made use of 



A , . . , .a ^ P 2 — P' 2 



<4 + w p-p' - <> ~ o,/ + ^p-p' ( 44 ) 



2M 

and 

|p 2 - p /2 I ^ |p - p'IIp + p' 



2M A ^ 2M A «IP-Pl<"P-P'> (45) 

analogously for and M B . 

Consequently, the non-relativistic solutions t/>(p) °f Eq. (^) are approximate solutions 
of the corresponding non-relativistic Schrodinger equation 



-f$ h+ (p - P->- v,(p,) = ( £ ' " Ma ~ Mb > {p) ' (46) 

or, written in the more familiar form in position space after a Fourier transformation, 

(-^-V 2 - a ^( r ) = (e' -M A - M b )^{v) . (47) 

As is well-known, in the case fi = the ground state solution of Eq. (EBf) is 



/ 32(aM r )5 1 

^ l(p) = V^^(^T^F (48) 

For a < 1, V'i(p) is strongly suppressed at p 2 > M 2 compared to p 2 <C Mj 2 , hence self- 
consistency is achieved. The wave functions of the excited states show similarly strong 
p-dependences. Note, however, that a priori other solutions of Eq. ([35]) with weaker p- 
dependences are not excluded for a -C 1, although they are certainly not expected on 
physical grounds. The numerical calculations presented in the next section therefore play an 
important role in the analysis of this limit. 

Another way to put the condition of strong p-dependence is to refer to the characteristic 
momentum scale aM r (from the expectation value (|p|)). Considering now the case /i ^ 0, 
it follows from phenomenological considerations (extension of the bound state) that the 
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characteristic momentum is given by the larger of aM r and /x. In particular, for the non- 
relativistic limit one needs that /i <C M r , in which case for the wave function a strong in- 
dependence qualitatively similar to Eq. is expected. As a result, the non-relativistic limit 
is self-consistent and is realized for a C 1 and /j <C M r , both conditions being necessary. 

The second limit of Eq. ([$5|) we will discuss, is the so-called one-body limit M B — > oo. In 
more physical terms, this means that M B 3> M\ , but also that M B is much larger than the 
characteristic momentum scale (p 2 ), which in turn is determined by the solutions ip(p) of Eq. 
(|35l) in this limit. The situation hence bears some similarity with the non-relativistic limit, 
and in fact one can use the analogues of Eqs. ( fP|) and ( f4"5| ) to determine the approximate 
form of the effective potential in this limit. As a result, in this situation Eq. ( p5|) can be 
approximated by 



M\ + p 2 



4>(p) + 



d 3 P ' 



V A (p,-p')i/>(p') = - M B )i>(p) 



(49) 



where 



V A (P,P') 



2M 



B 



2u)£2u$ 2uj p-p' 



uj. 



p-p 



p-p 



(50) 



Eq. ( f49D is a relativistic Schrodinger equation for particle A in the (static) potential Va- 

We will now show that Eq. (|4"9]) is identical with the equation for particle A interacting 
with a fixed source through the exchange of scalars of mass /i, thus establishing the inner 
consistency of the formalism in the one-body limit. To this end, consider the theory defined 
by the Hamiltonian 

H' = H' + Hj , 

H' = Jd 3 x 0^(x)(m 2 A -V 2 )0^(x) + ^(x)( / i 2 -V 2 Mx 



d 3 2#^(x)^(xMx) : + -f^(0) . 



2M, 



(51) 



B 



The form of Hj is motivated by considering 4>b{x) as a classical (positive-energy) Klein- 
Gordon field, with probability density 



d d 

P( x ) = <pB( X ) l -Q t <pB{x) ~ ^ B (x)l — (j) B (x) 



(52) 



In the case that M B is much larger than the relevant momenta p 2 , one has approximately 



p(x) = 2M B( f) B (x)(j) B (x) 
A particle localized at x = thus corresponds to 



(53) 



>b{ x )<Pb[x) 



2M, 



-5(x) . 



(54) 
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Figure 6: The second-order contributions in the one-particle sector for the interaction with 
a fixed source, represented by double lines in the diagrams. Only the last two diagrams lead 
to an interaction, given by Eq. (|55|). 

Using the latter formula in the classical counterpart of Hj in Eq. (||) leads to the interaction 
Hamiltonian Hj. 

Considering the one-particle states as in Eq. (|T9|), the Bloch- Wilson Hamiltonian cor- 
responding to H' leads, to lowest non-trivial order, to the diagrams represented in Fig. ^. 
The first four diagrams are diagonal in the basis {|pa)}- The first two, unlinked, diagrams 
reproduce the second-order correction E' n — E' to the vacuum energy,^ while the following 
linked but disconnected diagrams renormalize the mass of the particle to Ma as before. The 
last two diagrams in Fig. |] correspond to the interaction of the particle with the source, and 
translate to the expression 

(Pa\V a \ P ' a ) = ~i^J_ rfte-^y^A^O-^O-x^^x^^x) 

~*2Mb/ dte ~^ j ^^;(0,x)^(0,x)A F (0-t,x-0). (55) 
It is reassuring to see that this potential can be obtained by substituting 

<*(*. X )V£>*) = 2^*(*) (56) 

for the wave function of the 5-particle in the potential Eq. (|31| ) previously obtained for the 
two-particle sector. 

The corresponding Schrodinger equation is 

+ p^ V(P) + J |^<P|^|PW) =(E- i&Mp) . (57) 
Using the non-covariant form Eq. ( [TB|) of the Feynman propagators, it is readily shown that 

(p|Va|p') = ^4(p,P / ), (58) 

the potential Eq. (|50|) obtained in the one-body limit of Eq. (|35|) . Except for a constant 
shift in the energies, Eqs. ( |49| ) and ( pTj) are hence identical, thus establishing the consistency 
of the formalism in the one-body limit. 

4 Technically, the static source belongs to the vacuum, so that the second diagram in Fig. || which appears 
to be a self-energy correction of the source, in fact contributes to the vacuum energy. Diagrammatically, the 
double line is not counted as an external line, consequently the second diagram is considered unlinked. 
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4 Solutions 



In this section, we present the numerical solutions of the eigenvalue equation Eq. (|35|) . We 
have restricted attention to the case \x = (massless exchange particle), which represents 
both the physically most interesting situation and the most sensible test for the formalism. 
We also restricted the calculations to the s-wave spectrum. There is, however, no problem 
in considering higher angular-momentum states as well. Due to the spherical symmetry 
of the Bloch-Wilson Hamiltonian different angular-momentum states decouple. Since the 
potential depends, in terms of the angular variables, only on the angle 9 between the vectors 
p and p', one can employ a partial-wave decomposition 

oo 

V(p,p',cos0) = ^2ai(p,p')Pi(cos9) 

1=0 

l,m 

where p = |p|, p = p/p, and analogously for p'. Although rather cumbersome, analytical 
expressions for a; can be obtained for arbitrary values of I. Since the potential is diagonal in 
the angular momentum quantum numbers I and m, the eigenstates of the Hamiltonian will 
be eigenstates of the angular momentum. 

For eigenfunctions with zero angular momentum, the analytical integration of the angular 
variables in the potential term in Eq. (|35|) gives 

d 3 p' 1 1 v(Ip'I) 



- P < — P ' 

1 f°° dp' p'i){p') 



( 27r ) 3 j2u$2uB 2 ^p-p' w p + w p-p' 

+ P + P' ~ U v' 



log 



Up + \P - P'\ - ^v> 



(60) 



The integrand contains a logarithmic singularity, which, with the proper care, can be inte- 
grated easily. We use an integration method with a fine grid-size around the singularity. 
Furthermore, because of the wide separation of the relevant scales appearing in the prob- 
lem in the weak-coupling limit, we took special care to distribute the grid points over the 
full physical range. In order to determine the spectrum, we calculate matrix elements of 
the Hamiltonian with a number of smooth basis functions which go into a numerical eigen- 
value routine. In our calculations we have used two different bases which we specify in the 
following. 

The Coulomb basis is the best basis to analyze the eigenstates of the Bloch-Wilson 
Hamiltonian in the weak-coupling limit. The Coulomb basis consists of the eigenstates of 
the non-relativist ic Coulomb problem (n = 1, 2, 3, . . .), 



, , / 32n 2 1 ( {np/Kf-l . 

W)~\l nK s [{np/Kf + lf n - 1 \{np/KY + l ] ' 
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Ground-state wave function for a = 1 

Non-relativistic wave function for a = 1 

— ■ — Ground-state H wave function for a = 1 
■ - Non-relativistic H wave function for a = 1 



le-06 0.0001 0.01 1 100 

p[M r ] 

Figure 7: Double-logarithmic plot of the ground-state wave function for a — 1, compared 
to the non-relativistic results. The H wave function refers to a mass ratio Ma/Mb as in the 
hydrogen atom, the other wave function is for equal masses. 

where n = aM r , and is a Gegenbauer polynomial. Before calculating the spectrum 
in this basis, in order to improve convergence, we apply a similarity transformation to the 
effective Hamiltonian in Eq. (|35D , by means of which 

V(jp, p') — 1 V(p, p') J^p^, , (62) 

while the kinetic term remains unchanged. 

Away from the weak-coupling limit, the eigenstates are quite different from the corre- 
sponding Coulomb wave functions. Furthermore, the Coulomb basis is not a good basis to 
describe these states. The problem lies in the fact that all Coulomb wave functions have the 
same high-momentum tail, and only their low-momentum part varies. Beyond the weak- 
coupling regime, the high-momentum tail of the ground state of Eq. (|35|) turns out to be 
suppressed with respect to the Coulomb wave functions (see Figs. [?| and||). Adjusting the 
overall scale of the Coulomb wave functions cannot cure this problem. 

Even though in the Coulomb basis the energy converges rapidly with the number of 
basis states, one can show that it does actually not converge to an eigenvalue of Eq. (|35|) . 
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— Ground-state wave function for a = 2 

— Non-relativistic wave function for a = 2 
— Ground-state H wave function for a = 2 

■ - Non-relativistic H wave function for a = 2 
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Figure 8: Plot of the ground-state wave function as in Fig. [7|, for a = 2. 



Therefore, we have used a different basis, with a dependence on the two parameters (3 and 



K. 



0n(p) 



y/3(2w + l)(2/3- 1 

K 3/2 [(p/ K )3 + \f 



P 

- 1 n 



[(p/«) 3 + 1P" 



(63) 



P„ being a Legendre polynomial. This choice of basis leads to very satisfying results over a 
wide range of intermediate values of a. For most purposes the values k = aM r and (3 = 3/2 
turned out to be sufficient. The wave functions are normalized to unity, 



(64) 



We have tested our results by comparing H BW ip(p) to Eip{p) for the ground state, with 
the energy E and wave function ip(p) determined as described above. The deviations are 
typically below one percent, which demonstrates that our numerical result for both the energy 
and the wave function constitutes an excellent approximation. For a = 2, the corresponding 
deviation in the case of the Coulomb basis is of the order of 25%. We have used up to 40 
basis states, the changes in the deviation being minimal for more than 10 basis states. One 
can show that the difference {Hbw ~ E)ip(p) is orthogonal to the basis with an accuracy 
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Figure 9: The spectrum of binding energies Eb = E' — Ma — Mb (cf. Eq. (|35|) ) for s-states 
in the equal-mass case, compared to the ground state energies of the Wick-Cutkosky model 
H, the Hamiltonian eigenvalue equation in a Fock space truncation |)25||, the Regge theory 



predictions |26) , the light-front calculation [22 , and the Haag expansion results [24] in their 
domain of validity. 



of 10 -9 for our numerical integration, for both the Coulomb and the other basis. This also 
demonstrates that the problem with the Coulomb basis has nothing to do with the numerics. 

The spectrum, as in most relativistic approaches, shows weaker binding than is expected 
from the extrapolation of the non-relativistic formula (see Figs. [| and |10J) . Consequently, 
the wave functions peak at a lower momentum. More interestingly, the power-law behavior of 
the high-momentum tail of the wave functions changes. Whilst the Coulomb wave functions 
behave like pip(p) ~ p~ 3 , the wave functions found here beyond the weak-coupling regime 
fall off more rapidly, with pip(p) ~ p~ 3 ' 5 . The slope of the long-range tail in configuration 
space is unaltered by the relativistic corrections (see Figs. [7| and ||). 

We have also made an attempt to visualize the non-local potential Eq. ([57]). Given 
a solution ^(p) of the Schrodinger equation Eq. (|35|) , one can define an "effective" local 
potential Vx(r) in position space by 



V l (t)iP(t) 



d 3 p' 



np.p'Mp') 



J F.T. 



(e 1 - ^Jm\ + P 2 - y/NQ + jA Hp) 



(65) 



F.T. 



where F.T. denotes the Fourier transform to position space. The definition Eq. ( |65| ) makes 
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a 

Figure 10: The spectrum of binding energies for s-states as in Fig. |9|, for the hydrogen- 
masses case. 

sense at least for the ground state, given that the latter has no nodes. Unfortunately, it 
is numerically very difficult to evaluate the expression Eq. over the whole range of r. 
However, we could determine that, somewhat unexpectedly, at the origin r = 0, Vl tends 
towards finite values, numerically Vl(0) = —15.0 for a = 1 and Vl(0) = —4.0 for a = 2. 



5 Discussion 

The Bloch- Wilson Hamiltonian has a number of practical advantages over other relativistic 
approaches. Most notable is the absence of the energy eigenvalue in the potential part of 
the two-body bound-state equation. This has important computational advantages over 
formulations that have an energy-dependent integral kernel. We did not have to make any 
approximation for this result; it is a natural consequence of our approach. In the present 
formulation, the Hamiltonian is not Hermitian. Consequently, the left eigenstates are not 
identical to the right eigenstates, even though all eigenvalues turn out to be real. As a 
consequence of the similarity transformation back to the if -invariant subspace, the bound- 
state equation is expressed solely in terms of the two-particle wave function; all the internal 
workings of the interaction have gone in the Bloch- Wilson Hamiltonian. 

To the present order in the perturbative expansion, all formally divergent contributions 
can be identified with proper on-shell Feynman diagrams, hence the renormalization pro- 
cedure used in covariant Lagrangian perturbation theory can be applied to the effective 
Hamiltonian. This may not be possible at higher orders, nor can it be expected, of course, 
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when a non-perturbative approach for the determination of the Bloch- Wilson Hamiltonian 
is chosen. In these cases, other and possibly not manifestly covariant renormalization pro- 
cedures will have to be used. However, none of this is necessary in the present case. In 
particular, to the order considered the self-energy corrections merely serve to renormalize 
the masses. 

Led by the results of the FSR approach |14, ^TJ, it is nowadays believed that the Bethe- 
Salpeter equation suffers, apart from other unphysical features, also from a large underbind- 
ing. Our results lie closer to the non-relativistic values than the Bethe-Salpeter results and 
most other relativistic approaches (see Figs. |9| and |TU]), i.e., we find indeed a much stronger 
binding. Recently a discussion has emerged as to how the instability of the Wick-Cutkosky 
Hamiltonian could affect the bound state results [p2|— [p4|. At the present level of approxi- 



mation, there seem to be no consequences of this unphysical feature of the model, however, 
how our approach can be extended to be sensitive to the instability is a topic for future 
research. 

Not only is the non-relativistic limit properly recovered, as we were able to demonstrate 
both analytically and numerically, but also the one-body limit (RUf, where one particle 



becomes infinitely heavily, is consistent; the Schrodinger equation reduces to the equation for 
one particle interacting with a fixed source. Therefore, unlike the Bethe-Salpeter equation, 
the present approach can be used to study heavy-light systems, such as the hydrogen atom. 

Finally, we comment on the invariance of the results under Lorentz boosts. In princi- 
ple, bound states in a moving frame, i.e., with total momentum different from zero, can 
be calculated by solving the effective Schrodinger equation Eq. ([^), and the results can be 
compared to the relativistic energy-momentum relation. At any rate, a perturbative Hamil- 
tonian approach which by its very nature is not manifestly covariant, cannot be expected to 
maintain dynamical invariances, like boost invariance, exactly. As in the description of any 
non-perturbative phenomenon, some symmetries of the underlying theory will be violated at 
any level of approximation. A bound state in a moving frame will then slightly differ from a 
bound state at rest. This difference will become smaller with increasing order in the coupling 
constant to which the Bloch-Wilson Hamiltonian is calculated. However, we emphasize that 
in the present approach the expansion parameter is a Lorentz-invariant quantity, and the 
necessary renormalization is carried through in a covariant way, so that the violation of boost 
invariance might be expected to be small. To what extent boost invariance is actually broken 
at the present order and how much of it can be restored by calculating the Hamiltonian to 
a higher order is a topic of current research. 
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